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ABSTRACT 

We use deep VLT optical and near-IR spectroscopy and deep Spitzer/IRAC imaging 
to examine the properties of two of the most luminous Lya emitters at z = 5.7. The 
continuum red-ward of the Lya line is clearly detected in both objects, thus facilitating 
a relatively accurate measurement (10-20% uncertainties) of the observed rest-frame 
equivalent widths, which are around 160 A for both objects. Through detailed mod- 
elling of the profile of the Lya line with a 3-D Monte-Carlo radiative transfer code, 
we estimate the intrinsic rest-frame equivalent width of Lya and find values that are 
around 300 A, which is at the upper end of the range allowed for very young, moder- 
ately metal-poor star-forming galaxies. However, the uncertainties are large and values 
as high as 700 A are permitted by the data. Both Lya emitters are detected at 3.6 fj,m 
in deep images taken with the Spitzer Space Telescope. We use these measurements, 
the measurement of the continuum red-ward of Lya and other photometry to con- 
strain the spectral energy distributions of these very luminous Lya emitters and to 
compare them to three similar Lya emitters from the literature. The contribution 
from nebular emission is included in our models: excluding it results in significantly 
higher masses. Four of the five Lya emitters have masses of the order of ^ 10^ Mq 
and fairly high specific star-formation rates (> 10-100 Gyr~^). While our two Lya 
emitters appear similar in terms of the observed Lya rest-frame equivalent width, they 
are quite distinct from each other in terms of age, mass and star formation history. 
Evidence for dust is found in all objects, and emission from nebular lines often make 
a dominant contribution to the rest frame 3.6 ^m flux. Rich in emission lines, these 
objects are prime targets for the next generation of extremely large telescopes, JWST 
and ALMA. 
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1 INTRODUCTION 



* Based on observations obtained at the European Southern Ob- 
servatory using the ESO Very Large Telescope on Cerro Paranal 
through ESO programs 275.A-5012, 076.A-0553 and 080.A-0237 
and on observations taken with the Spitzer Telescope through GO 
programs 50308 and 60059 
f E-mail: clidmanOaao. gov.au 



Characterising the properties of very high redshift star- 
forming galaxies is important for a number of studies. For 
example, it is now generally accepted that ionising photons 
from massive stars in these galaxies led to an important 
event in the history of the Universe: the reionisation of hy- 



drogen in the intergalactic medium (IGM) ( Robertson et al 
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2OIOI). With the advent of Wide Field Camera 3 (WFC3) 



on the Hubble Space Telescope (HST), we are making rapid 
progress in discovering significant numbers of these galax- 



ies ([Bunker et al.||2010| |Bouwens et al.||2010| |McLure et al 



2010); however, progress on characterising their properties 



is slower, largely because these galaxies are extremely faint. 

Detailed spectroscopic follow-up of most of these ob- 
jects will need to wait for the upcoming generation of space 
and ground-based facilities, as the signal-to-noise ratios that 
are obtainable with current facilities are generally too low 



(e.g. Capak et al. 2011 Lehnert et al. 2010 Stark et al 



2010). Alternatively, one can study nearer and brighter ob- 



jects that are expected to be low-redshift analogues of these 
very high redshift sources (e.g. |Izotov et al.|2009 l, or, as we 
do in this paper, study the brightest high-redshift examples. 

Very high redshift galaxies (defined here as galaxies 
with redshifts greater than 5) are expected, intrinsically at 
least, to be powerful emitters in the Lya line (Hayes et al. 



2010). However, Lya is a resonance line with a large opac- 



ity, so even a small amount of neutral hydrogen modifies 
the line profile. Just how much and in what way depends on 
the intrinsic properties of the line and the distribution and 
kinematics of gas and dust in both the interstellar medium 
(ISM) attached to the galaxy and the IGM along the line-of- 
sight. Many high redshift galaxies do indeed emit strongly 
in the Lya line and are often called Lya emitters or LAEs 
for short. 

Between 2000 and 2002 we undertook a survey to find 
the brightest LAEs at z = 5.7 with the aim of furnish- 
ing a sample that could be studied in some detail with the 
new generation of telescopes that were becoming available at 
that time. This survey, called WFILAS (Wide Field Imager 
Lyman Alpha Survey), was conducted with the ESO 2.2m 
telescope and its Wide Field Imager (WFI), in conjunction 
with three narrowband filters and one intermediate band 
filter, all located in a region of low terrestrial background 
at ~ 815nm ( [Westra et "aI1|2006[ ). The survey covered 0.74 
square degrees and furnished 7 LAEs above a flux threshold 
of 5 X 10"^^erg/s/cm^ 

In this paper, we investigate the properties of two of the 
brightest LAEs from WFILAS survey using deep VLT op- 
tical and near-IR spectra and Spitzer/IRAC images. They 
are among the brightest LAEs currently known at redshift 
z = 5.7, and are therefore amenable to detailed study with 
current facilities. In Section 2 of the paper, we describe the 
VLT and Spitzer observations. In Section 3, we analyse the 
spectroscopic data, fitting the profile of the Lya line, mea- 
suring the continuum and searching for evidence of other 
emission lines. In Sections 4, we combine the spectra with 
the IRAC photometry to model the spectral energy distri- 
butions. We also reanalyse three other bright LAEs from 'Lai' 



et al. (20071. In Sections 6, we summarise out main findings. 



Throughout the paper, all equivalent widths are reported in 
the rest frame, unless explicitly noted otherwise, and we as- 
sume a fiat ACDM cosmology with Qa = 0.73 and Hq = 71. 



z = 5.7. From earlier observations (Westra et al. 2005 2006 ), 
the 2a lower limit on the EW was above 100 A for both ob- 
jects, which made them potential Population III galaxies and 
therefore interesting targets to search for Population HI sig- 
natures (cf. l Schaerer|2002 I. They are sufficiently bright that 
they can be observed in some detail with current instrumen- 
tation and before the availability of JWST, ALMA and the 
next generation of extremely large ground-based large tele- 
scopes. The properties of these LAEs, taken from [Westra et| 
al] p005l [2056^, are summarised in Table [l] 

Our observations consisted of deep medium-resolution 
spectroscopy with F0RS2 and SINFONI on the VLT 
and 3.6 and 5.8 ^m imaging with IRAC on the Spitzer 
Space Telescope. The FORS2 spectra cover the rest frame 
1180 to 1400 A region, which includes Lyman alpha and 
NVA1240. The SINFONI spectra cover the rest frame 
1640 to 2000 A region, which includes several lines, such 
as the CHI] AA1907,1909 doublet, and, for Sll 5236 only, 
HeIIA1640. At z = 5.7, the [OIII] AA 4959,5007 doublet and 
H/3 lie within the passband of the IRAC 3.6 /im filter. 



2.1 FORS2 spectroscopy 

SGP 8884 and Sll 5236 were observed with the multi-object 



spectroscopic (MOS) mode of F0RS2 (Appenzeller et al 
nigsl on Antu (VLT-UTl). The MOS mode of FORS2 con- 
sists of 19 movable slits that can be used to select targets 
over a 7' field. The slit length varies between 20" and 22". All 
data were taken with the FORS2 1028z holographic grating 
and the OG590 order sorting filter, which has a throughput 
of almost 80% in the region surrounding 810 nm. The slit 
width was set to 0"8. The details of the FORS2 observa- 
tions are listed in Table [2| 

One MOS slit was placed on the LAE. Other slits were 
placed on a selection of galaxies and stars. The stars were 
used to monitor the spectral PSF. Individual exposures were 
of 1370 seconds duration, and two of these were taken in a 
single observing block (OB). Several of these OBs were exe- 
cuted over the course of several nights to reach the necessary 
depth. To facilitate the removal of bright sky lines and to 
reduce the impact of cosmic rays and CCD charge traps, the 
object was placed in a slightly different part of the slit for 
each exposure. 

The data were reduced in a standard manner using 
irafQ tasks. The bias was removed by fitting a low order 
polynomial to the overscan region and pixel-to-pixel varia- 
tions were removed with lamp fiats. The removal of night sky 
lines, cosmic rays and the co-addition of the 2-dimensional 
spectra were done using our own software. These sky sub- 
tracted and co-added 2-dimensional spectra are shown in 
Fig. [1] Although the continua are considerably fainter than 
the Lya emission, they are still visible in the 2-dimensional 
spectra (Fig.[TJ. The 1-dimensional spectra of the LAEs were 
then extracted and calibrated in wavelength and fiux. 

The signal-to-noise ratio of the Lya line varies from 100 
for Sll 5236 to 150 SGP 8884. The integrated Lya fiuxes are 



2 TARGET SELECTION AND OBSERVATIONS 

With luminosities of ~ 3 x 10'^^erg s~\ SGP 8884 and 
Sll 5236 fi-om the WFILAS survey of [Westra et~aL] ( [2006] ) 
are among the most luminous LAEs currently known at 



^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories which are operated by the Association of Universities 
for Research in Astronomy, Inc., under the cooperative agreement 
with the National Science Foundation 
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Table 1. Properties of SGP 



and Sll 5236 



Object 


RA 


Declination 




Flux (imaging)^ 


Flux (spectroscopy)^ 








[A] 


[10"-'-'' ergs"-*- cm~^] 


[10"-'-'' erg s"-' cm~^] 


SGP 8884 


00:45:25.38 


-29:24:02.8 


> 125 


9.5 ± 1.4 


9.05 ±0.06 


Sll 5236 


14:43:34.98 


-01:44:33.7 


> 100 


7.0 ± 1.2 


4.18 ±0.04 



Westra et al. 



VVestra et al. 



20061 



Note 1: Initial 2(t lower limits from 
Note 2: Narrow band fluxes from ' 
Note 3: Integrated line flux from this paper. Does not include slit losses 



1 (1200 5 
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Figure 1. Sky subtracted spectra of SGP 8884 (top two panels) and Sll 5236 (bottom two panels). The spectra are split into two for 
display purposes. Noisy regions mark the location of bright night sky lines. The continuum red-ward of the Lyman alpha line is visible 
in both objects, although it is considerably fainter for Sll 5236. The spectral trace of an [Oil] emitting galaxy lies about 6" below the 
spectral trace of SGP 8884. At a resolution of 5000, the [Oil] AA3727 doublet is clearly resolved. 



Table 2. FORS2 Observations 



Object 


Total exposure time 


Median seeing 
["] 


SGP 8884 
Sll 5236 


16440 
43840 


0.65 
0.77 



within a factor of two of the fluxes derived from narrow-band 
imaging ( |Westra et al.|2006[ ). The details are listed in Table 
[l] The larger difference for Sll 5236 can be attributed to the 
poorer observing conditions during which the spectra were 
taken. In the remainder of this paper, we scale the FORS2 
spectra so that the Lya fluxes computed from spectroscopy 
match those derived from narrow-band imaging. 



2.2 SINFONI spectroscopy 

SGP 8884 and Sll 5236 were observed with SINFONI on 
Yepun (VLT-U T4). SINFONI | |Eisenhauer et al.|2003||Boir 
|net et al.|2004l ) is a near-infrared integral field spectrograph 
that has an adaptive optics (AO) module that can be used 



with natural guide stars or a laser guide star to improve the 
spatial resolution. The observations were done without AO 
due to absence of sufficiently close and bright field stars for 
AO use. 

The targets were observed with the 0'.'25 objective, 
which has a field of view of 8" x 8" and 0'.' 125 x 0'.'25 spaxels, 
and the J grism, which covers the 1095 nm to 1350 nm wave- 
length range with a resolution of approximately 2000. Each 
observing block consisted of four 600 second exposures. Be- 
tween exposures, the telescope was offset by a small amount 
so that the data could be used to estimate the sky while the 
object remained within the SINFONI 8" field of view. The 
details of the SINFONI observations are listed in Table |3l 

The data were reduced following the steps described 



in Cuby et al. (20071. The spectra were extracted with 1 
apertures that were placed at the locations the objects were 
expected to be. The spectra and the associated error spectra 
were calibrated in wavelength and flux. 

The continuum of neither object was detected, even af- 
ter binning the extracted Id spectra. No emission lines were 
detected either. However, the SINFONI data allow us to 
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Table 3. SINFONI Observations 



SGP 8884 » IF^C 



Object 


Total exposure time 


Median Seeing 

["] 


SGP 8884 
Sll 5236 


24000 
9600 


0.8 

0.7 



place upper limits on the most likely emission lines that lie 
in the region covered by the SINFONI data (see Sec. |3.3[). 



2.3 IRAC imaging and photometry 

SGP 8884 was imaged at 3.6 and 5.8 ^im with IRAC on the 
Spitzer Telescope during cycle 5, a period during which the 
instruments were still being cryogenically cooled. To reach 
the required sensitivity, the total integration time was set to 
36,000 seconds, which was split into 2 AORfj^ Each AOR 
consisted of 180 100-second exposures. Given the large num- 
ber of offsets and the characteristics of the target, we chose 
the cycling dither pattern with a medium scale factor. 

Sll 5236 was also imaged with IRAC. However, the ob- 
servations were carried out during cycle 6, when Spitzer had 
run out of cryogens, so only the 3.6 /im data could be taken. 
The integration time and dither pattern were identical to 
that used for SGP 8884. 

Individual integrations were processed using versions 
S18.7.0 and S18.18.0 of the basic calibration data pipeline 
for SGP 8884 and Sll 5236, respectively. The processed data 
were then combined to an image with a final pixel scale of 
0'.'61, a factor of two finer than the real pixel scale of the 
IRAC cameras, using version SI 8. 3.1 of the MOPEX mosaic- 
ing software. 

For SGP 8884, the IRAC 3.6 and 5.8 ^tm images reach 
1 (7 point-source sensitivities of 0.03 /^Jy and 0.32 /ijy at the 
location of the target. The 3.6 ^m image is shown in Fig. [2] 
together with a F0RS2 image in the 815-13 narrow band 



filter that was used to define the MOS mask (Sec. 2.1 1. Both 
images are centred on SGP 8884. SGP 8884 is undetected 
in the IRAC 5.8 /xm image. 

For Sll 5236, shown in Fig.[3j the 1 a point-source sen- 
sitivity of the 3.6 /im image is 0.08 /iJy. Since Sll 5236 was 
observed during the warm mission, there were no 5.8 /^m 
data. 

Both SGP 8884 and Sll 5236 were selected because they 
were very bright LAEs and because they lie in relatively un- 
crowded regions. This is quite important, as the sensitivity 
of the IRAC 3.6 micron observations is well below the limit 
at which confusion from overlapping sources becomes a seri- 
ous issue. The centroids of SGP 8884 in the IRAC [3.6] and 
FORS2 815-13 images are consistent with each other, sug- 
gesting that there is little contamination from other objects 
that happen to land on the line-of-sight. It also means that, 
to within the uncertainties in measuring centroids (~ 0.5 kpc 
at this redshift), the regions responsible for the Lya emis- 
sion and the 3.6 flux are coincident. On the other hand, for 
Sll 5236, there is an offset of approximately 0'.'5 (~ 3 kpc) 



Astronomical Observing Request 




F0RS2 

815/13 nm 



10' 



Figure 2. IRAC 3.6^111 (left) and FORS2 narrow band (right) 
images of SGP 8884. The images are 45"across with North is up 
and East is left. The FORS2 image was taken with the 815-13 
narrow band filter and was used to prepare the FORS2 MOS 
masks. Most of the flux in the FORS2 image is from Lya. 



Sll 5236 IRAC V 

W 3.6 micronE 



•1 



F0RS2 

815/13 nm 



iO' 



Figure 3. As for Fig.|2] but for Sll 5236. Most of the flux in the 
FORS2 image is from Lya. 



between the two centroids. This may mean that there is a 
faint object, not seen in the F0RS2 image, biasing the cen- 
troid in the IRAC 3.6 ^m image or that the region that is 
responsible for producing the flux at 3.6 /im is offset from 
the from the region that produces the Lya emission. 

Fluxes were measured in 2"radius apertures. The centre 
of the aperture was determined by the position of LAEs in 
the F0RS2 images relative to sources that were common to 
and clearly detected in both the FORS2 and IRAC images. 
The fluxes were then corrected to a larger aperture (radius 
4") using an aperture correction that was determined from 
relatively bright stars. The fluxes are reported in Table |4] 

The flux uncertainties are computed from the uncer- 
tainty images produced by MOPEX, which includes only sta- 
tistical errors. Systematic uncertainties, such as the error 
that comes from unresolved sources, are not included. 



Table 4. Fluxes in the observer frame 



Object 


0.87 /im flux 


3.6 /im flux 


5.8 /im flux 






(MJy) 


(MJy) 


SGP 8884 


0.23 ± 0.04 


0.52 ± 0.06 


0.32" 


Sll 5236 


0.17 ± 0.04 


0.48 ±0.12 





"■la upper limit 
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2.4 Morphologies 



3.2 Lya fitting 



Neither SGP 8884 nor Sll 5236 are clearly resolved in tiie 
FORS2 images tliat were used to prepare tlie maslcs for spec- 
troscopy. The image quality was around O'.'S in both images, 
so the Lya emitting regions are smaller than 3 kpc. 

However, a field containing Sll 5236 was observed with 
the Advanced Camera for Surveys instrument (ACS) on 
HST as part of program 10798 (PI Koopmans) with both 
the F555W filter, where it is not detected, and the F814W 
filter, where it is clearly is detected. The FWHM of a Gaus- 
sian fit to the profile of Sll 5236 is 0'.'2, which is consider- 
ably broader than the PSF and corresponds to physical size 
1.2 kpc at the redshift of Sll 5236. 

Approximately half of the fiux in F814W ACS image of 
Sll 5236 comes from the Lya line. The rest comes from the 
continuum red-ward of the line. 



3 ANALYSIS OF THE SPECTRAL 
PROPERTIES 

3.1 Modelling the continuum and deriving the 
rest-frame EW of Lya 

By inspecting the sky subtracted spectra in Fig.[T] it is clear 
that the continuum red-ward of the Lyman alpha line is de- 
tected in both objects. We fit a power law with index /3 to 
the data red-ward of the line and use the fit to estimate the 
fiux of the continuum at Lya and the observed rest-frame 
equivalent width (EW) of the line. The wavelength inter- 
val of the fit is indicated in column 2 of Table [5] A power 
law provides an adequate fit to the data, although given 
the uncertainties, the error in /3 is substantial and allows 
one to constrain neither the age of the stellar population re- 
sponsible for the UV emission nor the amount of reddening. 
Excluding the SINFONI data from the fit does not change 
the value of /3 significantly. 

Blue- ward of Lya, the continuum is clearly suppressed 
relative to the continuum red-ward of the line. In fact, we 
do not detect a significant signal in either object. We fit the 
average level of the continuum on the blue side of the line 
using all the data within the range indicated in column 6 of 
Table m The mean values for SGP 8884 and Sll 5236 are 
consistent with zero within 2-a. 

The ratio between the measured continuum between 
Lya and Ly/3 and the one extrapolated from the continuum 
red-ward of Lya is often used to measure the fraction of neu- 
tral hydrogen in the IGM. The ratio has been well measured 
for QSOs (Songaila 2 004) , and at ^ = 5.7 is around a few 
percent, which corresponds to a neutral fraction of ~ 10"*, 



with some scatter between objects (Fan, Carilli & Keating 



2006) 



The upper limits for our two LAEs are a factor of five 
less constraining. If we extrapolate the continuum red-ward 
of Lya to 1185 A (the upper end of the wavelength range 



used in Songaila (2004l), we find a flux ratio of 0.23 ± 0.11 
for SGP 8884 and 0.14 ± 0.11 for Sll 5236. While this in- 
dicates that the continuum blue-ward of the line is affected 
by neutral hydrogen in the IGM, it does not show that the 
line itself is affected. We will return to this issue later in this 
section. 



The spectra of both LAEs, centred on Lya, are shown in 
Fig. |4] Both lines are clearly asymmetric, with steep in- 
creases on the blue side of the lines and long red tails. Very 
likely, the asymmetry is indicative of an expanding ISM 



(Verhamme, Schaerer, & Maselli 20061. A relatively weak 
shoulder on the red side Lya line of Sll 5236, tentatively 



identified in Westra et al. (20051, is now clearly visible. We 



fit the line with a variety of models, first trying analytic 
models that have been used in the past and then more so- 
phisticated 3-D Monte Carlo radiative transfer models. 



3.2.1 Analytic models 



Following Dawson et al. (20021; Hu et al. (20041 and Westra 



et al. (2005 20061, we fit the profile of the Lya line with a 
truncated Gaussian that is convolved with the instrument 
resolution. There are four parameters in the fit: the FWHM 
and centre of the truncated Gaussian, the instrument resolu- 
tion, which is itself modelled as a Gaussian, and the overall 
normalisation. Since the average seeing during the observa- 
tions was narrower than the width of the slit, the instrument 
resolution is mostly defined by the seeing. The continuum 
blue-ward of the line centre is set to zero. Red-ward of the 
line, the continuum is set to the value determined in Sec. |3.1| 
The fits are shown in Fig. [4j and the best fit values for the 
parameters of the model are shown in Table [6] The fits are 
constrained using the entire wavelength shown in Fig. |4] 

On the blue side of the line, the models fit the data 
very well. The instrument resolution, which is determined 
during the fit, is around {R ~ 5000) and agrees with the 
reported seeing. The blue side of the line is consistent with 



a vertical cutoff. Ouchi et al. (20101 reported the detection 



of a "knee" on the blue side of the profile in LAEs at 2 = 5.7 
and z = 6.5. Our data have considerably better signal-to- 
noise and better resolution (~ 1.6 A). There is no evidence 
for such a feature in the two LAEs that we observed. 

On the red side, however, the models are a poor fit 
and fail to reproduce the level of detail seen in the data. 
For Sll 5236, the shoulder at ~ 8183 A is, not surprisingly, 
completely missed. The fit for SGP 8884 is better; however, 
there also appears to be excess fiux in the red wing of the 
line around ~ 8105 A. Such features are characteristic of an 
expanding dusty ISM (cf. Verhamme et al!||2008 l. We now 
move to using a 3-D Monte Carlo radiation transfer code to 
model escape of Lya photons from the LAEs and to see if 
we can reproduce the red wing better. 



3.2.2 3-D Lya radiation transfer 

To model the Lya profiles, we use an updated version of 
the 3-D Lya radiation transfer code MCLy a of [Verhamme, 



Schaerer, & Maselli ( 2006 \ described in Schaerer et al 



|2011 


I. The code has been used successfully in 


Verhamme 


et al. 


(20081; Schaerer & Verhamme (20081; 


Dessauges- 


Zavadsky et al. (20101; Vanzella et al. (20101 to model Ly- 



man break galaxies (LBGs) and LAEs. We follow [Verhamme| 
et al. ( 2008) and fit a simple model to the data. The model 



consists of a line and continuum emitting region that is sur- 
rounded by an expanding, dusty shell of neutral hydrogen. 
The dust is assumed to be uniformly distributed within the 
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Table 5. Quanitites derived from to the fits to the continuum red- ward and blue- ward of Lya. 



Object 


Fit interval 


/3 


Continuum flux at Ly« 


EW 


Fit interval Continuum flux below Lya 




[A] 




[le-20 ergs/s/cm^/ A] 


[A] 


[A] 


[le-20 ergs/s/cm^/ A] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


SGP 8884 


1222-1797 


-1.7± 1.5 


8.2 ± 1.0 


166 ±20 


1169-1214 


1.9 ±0.9 


Sll 5236 


1243-1778 


-3.7 ±2.9 


3.9 ± 0.8 


160 ± 36 


1147-1214 


0.6 ±0.5 




-0.2 t 



8080 



Observed Wavelength [A] 



8140 



8180 



Observed Wavelength [A] 



Figure 4. Model fits to SGP 8884 (left) and Sll 5236 (right). The analytic model used in the fits is shown as the blue solid line. 
VLT/FORS2 data are plotted as red points. 



shell. There are seven parameters in the model. In practice 
we use an automated fitting tool, which synthesises Lya pro- 
files from the extensive library of more than 6000 radiative 



transfer models of Schaerer et al. (20111 covering a wide 



range of parameter space. The continuum at Lya is fixed to 
the value listed in Table [5] as these have been determined 
over a larger wavelength region than the regions used to fit 
LyQ and is therefore more accurate. 

At z — 5.7, it is possible that resonant scattering of Lya 
photons by neutral hydrogen external to the galaxy - either 
in the form of nearby Lya forest systems or remnant neutral 
gas in the IGM - may affect the Lya profile. Stochastic as 
this is, it provides a Monte Carlo problem all of its own, 
and one that we do not have the possibility to address on 
a case-by-case basis. Instead, we adopt two simulations per 
galaxy, which are designed to bracket a broad range of pos- 
sibilities: we firstly assume that the IGM does not affect the 
lin^ and secondly assume that the IGM absorbs all the 
fiux blue-ward of the line centre. In this second approach 
the intervention occurs after radiation transport (to main- 
tain the effects of the ISM) but before convolution with the 



This does not exclude the possibility that the IGM starts to 
suppress the continuum at some point that is blue-ward of the 
line. 



instrumental profile, and thus is somewhat comparable to 
the method used in Sec. |3.2.1[ but with the shape of the 
red side determined by simulation. We present the results 
from the first simulation first. The best fit parameters and 
the resulting escape fractions of Lya photons and the qual- 
ity of the fits as characterised by the reduced &re shown 
m Table ID The escape fraction is defined as the fraction of 
Lya photons that escape the galaxy. It does not include the 
fraction that are scattered by the IGM. 

The redshifts are well constrained and we use these val- 
ues throughout the remainder of the paper. Column densi- 
ties of neutral hydrogen are determined to be small, just 
0.3 X 10^® cm-2 and 2 x 10^^ cm"^ for Sll 5236 and 
SGP 8884, respectively. These values are unconstrained at 
the low Af(HI) end, but at the high end are strongly con- 
strained, with statistical upper error-bars of 0.4 dex. The 
characteristic features that allow more accurate determina- 
tion of effective dust optical depth, t, do not appear in ei- 
ther of these spectra, so the amount of absorption is poorly 
constrained, although reasonable lower-limits are obtained 
for both galaxies: we determine lower limits of r = 0.5 for 
SGP 8884 and 0.7 for Sll 5236 at the 68.3% confidence 
limit, although the upper limit is unconstrained by our grid, 
even at the 68.3% confidence limit. For both galaxies the 
full range is permitted at the 99.7% confidence limit. 

The Lya escape fractions are found to be small for 
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Table 6. Line fitting results. Columns 3 through to 9 are obtained from the fits and are, respectively, the redshift (z), the expansion 
velocity of the HI shell (\4xp), the column density of neutral hydrogen in the shell {N{H)), the amount of dust absorption (r), a Doppler 
parameter describing the random motions of the neutral gas (6), the intrinsic FWHM of the Lya line (FWHMq), and the intrinsic 
equivalent width of the Lya line (EWq). The escape fraction /esc, is derived from the fitted parameters. 



Object 


Fit Type 


z 


Vexp 


logAf(HI) 


T 


b 


FWHMo 


EWo 


/esc 










[kms-i] 


[cm^^] 




[kms^^] 


A 


A 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


SGP 8884 


Analytic. 


5.6527 
















2.1 


SGP 8884 


MCLya 


5.6517 


250 


19.3 


3.0 


80 


300 


350 


0.0232 


2.71 


SGP 8884+IGM 


MCLya 


5.6517 


50 


16.0 


2.0 


40 


600 


350 


0.104 


1.68 


Sll 5236 


Analytic 


5.7242 
















2.9 


Sll 5236 


MCLya 


5.7245 


300 


18.5 


1.0 


80 


200 


250 


0.300 


2.03 


Sll 5236+IGM 


MCLya 


5.7237 


400 


18.0 


0.2 


80 


300 


350 


0.819 


1.01 



SGP 8884 and moderate for Sll 5236, taking values of 2.3 
and 30%, respectively. For SGP 8884 this value is not tightly 
constrained and /esc = 50% (100%) are permitted at the 
68.3% (99.7%) confidence Umit. However for Sll 5236 /osc 
is constrained to < 50% at the 99.7% confidence level. 

The intrinsic rest-frame equivalent widths of the Lya 
line is 350 A for SGP 8884 and 250 A for Sll 5236, with 
values as high as 400 and 500 A permitted for Sll 5236 and 
SGP 8884, respectively (68.3% confidence). The lower limits 
just reproduce the measured EWs given in Table [5] 

The radiation transfer code (the blue curve in Fig. |5| 
provides a better fit to the the main features of the Lya 
line in Sll 5236 than the analytic model used in Sec. |3.2.T| 
In particular, the red shoulder that is seen in Sll 5236 is 
reproduced by the model. Conversely, for SGP 8884, the 
radiation transfer code provides a poorer fit before the IGM 
is considered. 

For both LAEs, the simulations suggest that at wave- 
lengths short-ward of the Lya line (~ -700 to -800 kms~^) 
the effects of radiation transport on continuum photons be- 
come negligible, and thus (modulo the dust extinction law) 
should behave no differently from photons on the red side 
of Lya (see Fig. [5|. At this scale, i.e. close to the Lya line 
the model is compatible with the observations. 

The second of our simulations assumes that all the 
transmitted flux blue-ward of the systemic velocity is ab- 
sorbed by the IGM. Obviously this removes all flux, includ- 
ing the small bumps, on the blue side of the line for both 
galaxies (magenta line cf. blue line). However, even with the 
quality of the present spectra it is not possible to distinguish 
from the observations blue-ward of Lya whether or not the 
IGM affects directly the Lya line. Both models fit the blue 
side of the Lya line equally well. In principle, the IGM may 
also affect the main red peak by giving the impression that 
the peak itself has shifted red- ward. This effect, anticipated 



by, for example, Haiman (20021; Santos (20041 and Dijk- 



stra fc Wyithe (20101 and noted in Verhamme et al. (20081, 



appears to be largely mitigated by the fact that the line is 
already shifted to the red from the multiple scatterings that 
occur in the ISM. 

Quantitatively the quality of the overall line fit is im- 
proved for both objects by the addition of the IGM, bringin; 
per degree-of-freedom down to 1.7 and 1.0 for SGP 8. 



^ For SGP 8884, there is some evidence from the scatter of the 



and Sll 5326 respectively and more in line with the values 
that we are looking for in this kind of analysis. The treat- 
ment of the IGM has further noticeable effects on the de- 
rived parameters that are worthy of discussion. Firstly it 
substantially reduces the dust optical depth of the ISM in 
both galaxies, although the formal level of constraint is no 
better than it was in the case where the IGM was ignored. 
In the ISM only model, the only way to remove emitted 
photons from the system is to absorb them by dust. Adding 
the IGM introduces an additional sink of photons. Thus, 
when we introduce the step-function IGM, we see higher es- 
cape fractions and consequently lower measures of r. Indeed 
the escape fractions have increased substantially, from 2.3 
to 10% in SGP 8884 and from 30 to 80% in Sll 5326. In- 
terestingly, the best fit intrinsic rest frame EWs are largely 
unchanged. The 68.3% confidence interval covers a slightly 
larger range, from 150 to 700 A. 

The column density of neutral hydrogen, A'^(HI), is sub- 
stantially reduced for both LAEs; however, the 68.3% con- 
fidence interval brackets the best fit values of the first sim- 
ulation. In the case of SGP 8884 the intrinsic FWHM of 
Lya doubles to 600 km/s; however, the uncertainty in this 
parameter is large. The 68.3% confidence interval extends 
down to 250 km/s. Such broad lines have been reported 



in other LAEs (e.g. Verhamme et al. 2008 Vanzella et al 



20101. They are more commonly observed in AGN; however. 



we see no evidence for AGN activity in either SGP 8884 or 
Sll 5236. The NVA1240 line, which is commonly seen in 
AGN, was not detected in either LAE. Upper limits on NV 
are listed in Table [T] 



3.2.3 Discussion 

Neither of the IGM models used in this paper are likely to 
represent the average state of the IGM along the line-of-sight 
to LAEs at z = 5.7. The first model assumes no attenuation 
of Lya from the IGM. The probability that this occurs for 



both objects is quite small (Laursen, Sommer-Larsen, & Ra- 
zoumov|2011 1. It is worth noting; however, that both LAEs 
were selected as being the brightest LAEs over a relatively 
large area of sky, so there will be a bias in selecting LAEs 



data in the continuum regions that the errors have been underes- 
timated by 20-30%. 
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that have more transparent hnes of sight than the average. 
The second model assumes that all photons blue-ward of 
the systemic redshift are scattered by the IGM. This step 
function neglects the damping wing of Lya, so the attenu- 
ation may extend to longer wavelengths depending on the 
kinematics of the IGM with respect to the source and the 
ionization state of the circum-galactic medium (see, for ex- 
ample,|H aiman'20 02[|Santos|2004l|Laursen, Sommer-Larsen,| 
fc Razoum ov 201]J7^ 

A more precise treatment of the IGM (as in, for exam- 
ple, Laursen, Sommer-Larsen, & Razoumov 2011 1 is not war- 



Table 7. Line ratios of other lines with respect to the observed 
Lya flux. Upper limits are 3-cr. 



ranted with the current data. One clear effect from adopting 
more accurate IGM models is that some of the fitted pa- 
rameters would almost certainly change. For example, if the 
IGM were to absorb flux red-ward of the systemic velocity 



(as it does in the models of Laursen, Sommer-Larsen, & 
Razoumov||2011 1, then the best fit redshift would likely de- 
crease, as it did for Sll 5236. Hence, an direct measure of the 
redshift from other spectral features, such as nebular lines 
in the near-IR, would be useful, as this would reduce the 
number of fitted parameters. Unfortunately, our attempts 
at detecting such lines with SINFONI proved unsuccessful 
(Sec 



3.31 



Similarly, a measurement of the expansion velocity of 
the shell from the velocity offset of low-ionisation interstellar 
lines would reduce the number of parameters further. Our 
data were not deep enough to see these lines. 

3.3 Other emission lines 

In addition to Lya, the F0RS2 and SINFONI spectra 
cover a number of lines commonly observed in AGNs and 
LBGs. These include, NV A1240, GUI] AA1907, 1909 and, for 
Sll 5236 only, HeIIA1640. The spectra do not cover the 
NIV] AA1483, 1487 doublet, which has been seen in some 
star-forming galaxies (iVanzella et al. 20101. Unlike Lya, 



these lines are relatively unafi'ected by the ISM and IGM, 
and can be used to derive redshifts. 

Apart from Lya, we do not detect any line emission in 
either SGP 8884 or Sll 5236. To derive upper hmits to the 
fluxes and to the line ratios with respect to the observed Lya 
flux, we compute the noise in 9 A bins centred at the red- 
shifted wavelengths of these lines using the redshifts derived 
from the fltting of the Lya line in Section [3.2| In Table [7] we 
quote 3-cr upper limits. The sensitivity limits vary strongly 
between lines because some lines overlap regions of bright 
telluric emission from OH and O2, while others are in rela- 
tively clear regions. 

In addition to listing GUI] and Hell, we provide 3-a 
upper limits for OIII] a nd NV. CHI], OIII] and Hell were 
detected in the Lynx arc ( Fosbury et al.|2003 1 and in Q2323- 
BX418 ( |Erb et al.|[2010| , a young, low metallicity, unred- 
dened galaxy at z=2.3, while NV is commonly found in 
AGNs (Vanden B erk et al.||2001[ ). While both the Ly nx arc 
and Q2323-BX418 are LAEs, they appear to be quite differ- 
ent to the average LAE. In addition to being rich in emission 
lines, the ionisation parameter for both objects is quite high, 
U = -1. 

In TablejT] we list the flux ratios between Lya and Hell, 
GUI] and OIII] in the Lynx arc. If the properties of SGP 8884 
and Sll 5236 were similar to the Lynx arc, then both GUI] 
and OIII] would have been clearly detected in SGP 



Line 


-^rcst 


bli OZJD 


DljF 8884 


Lynx 


NV 


1240 


< 0.043 


< 0.010 


< 0.084 


Hell 


1640 


< 0.198 




0.033 


OIII] 


1661 


< 0.191 


< 0.109 


0.058 


OIII] 


1666 


< 0.133 


< 0.042 


0.121 


GUI] 


1907 


< 0.148 


< 0.085 


0.113 


GUI] 


1909 


< 0.124 


< 0.106 


0.077 



and marginally detected in Sll 5236. In this respect, both 
SGP 8884 and Sll 5236 are different to the Lynx arc and 
more similar to the average LAE. 

In Sec. [4] we flt SED models that include nebular emis- 
sion lines to the broadband photometry. We can check to 
see if the upper limits listed for the lines in Table [7] are con- 
sistent with the strength of these lines in best flt SEDs. For 
both Sll 5236, where the lines in the SED model are very 
weak, and for SGP 8884, where the lines are stronger, this 
is indeed the case. 



4 SED FITTING 

To gain further insight into the properties of our objects 
and other LAEs, we have carried out SED fits using the Hy- 



per z code (Bolzonella, Miralles, & Pello 20001 modified by 



Schaerer & de Barros ( 2009 1 to take into account the effects 
of nebular emission (lines and continua) on broad band pho- 
tometry. These effects are expected to be largest for objects 
with emission lines (as the case for the LAE discussed here) 
at high redshift, since the contribution of emission lines to 
photometric bands increases with (1 + z). The code has al- 
ready been successfully applied to several samples (see e.g. 
Schaerer fc de Barros|2010l|Dessauges-Zavadsky et al.|2011[ ), 
and the method has also been implemented by other groups 
and applied to different samples, in particular to LAEs at 
z ~ 2-3 and at 6 ( |Ono et aL]|2010[ [Acquaviva, Gawiser, fc| 
Guaita|20l"T| . 

For the SED fits we use the spectral templates of 
[Bruzual fc Gharlot| (|2003[) to which we add nebular emission 
as described in Schaerer & de Barros (20091. We consider 
three metallicities between solar and 1/200 Zq, exponen- 
tially decreasing star-formation histories (oc exp"*''^) with 
r varying between lOMyr and 00 (constant star- formation) , 
and variable extinction described by the Calzetti law. We 
have also considered models with a fixed, constant tar- 
formation rate (SFR). Non-detections are included in the 
minimisation by setting the observed flux to zero and tak- 
ing the la upper limit as the error. The SED fits gives the 
resulting fit parameters, age, stellar extinction -Eb-v , stel- 
lar mass Mi, SFR. The probability distribution function of 
these parameters (and derived uncertainties) are determined 
by Monte Garlo simulations, generating a large number (typ- 
ically 1000 or more) of realisations of the SED. The stellar 
masses and SFRs quoted here are based on a Salpeter IMF 
from 0.1 to 150 Mq. For comparison with a more realistic 
IMF including the observed turn-over at low stellar masses 
these values must be corrected downward (e.g. by a factor 
1.8 for a Chabrier-type IMF). 
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Figure 5. Model fits to SGP 8884 (left) and Sll 5236 (right). The model in which the IGM attenuation is ignored is shown by the 
blue unbroken line, while the magenta line shows the model for which the transmission through the IGM was set to zero blue-ward of 
the systemic redshift. The intrinsic spectra (i.e. before radiation transport and attenuation) is shown by the dotted line in both cases. 
VLT/FORS2 data are plotted as red points. The vertical dashed lines mark the line centre of Lya for the four cases. For SGP 8884, the 
line centre of Lya overlaps in the two cases considered. The velocity scale in the right hand figure is only valid for the case where IGM 
attenuation is ignored. The insets illustrate magnified views of the continuum blue-ward of Lya, approximately within 1000 km s^^ of 
line centre. 



4.1 SED modelling of SGP 8884 and Sll 5236 

The broadband and spectroscopic measurements of the 
SED, summarised in Tabl e [4| and the upper limits for the 
B and R bands from [Westra et al.| ( |2006| are used to con- 
strain the SED fits. Given the relatively small number of 
observational constraints it is not surprising that good SED 
fits are easily obtained (see Figs. [6] and [7|. We do not use 
the strength of Lya that is derived from the modelling of the 
line to constrain the fits directly, as the uncertainty in the in- 
trinsic strength of Lya is too large. For Sll 5236, however, 
most models with variable star-formation histories do not 
predict a strong enough Lya emission. Therefore we have 
also examined fits imposing a constant SFR. 

The main qualitative difference between the two best-fit 
solutions obtained with templates including nebular emis- 
sion is obviously the importance of the emission lines in 
the rest-frame optical, which is predicted to be stronger for 
SGP 8884 than for Sll 5236. The reason is the non-detection 
at 5.8 /xm in SGP 8884 with a fiux below that at 3.6 /im, 
which drives the solution to a younger age, hence stronger 
emission lines. For SGP 8884, approximately two-thirds of 
the 3.6 fim fiux comes from line emission. 

The probability distribution function (pdf) for the 
age, attenuation, stellar mass, SFR, and specific SFR {— 
SFR/ Mi,) for the two objects is shown in Fig. [o] (black and 
red lines). The effect of imposing a constant SFR is shown 
by comparing the solid and dashed lines. As can be seen 
from these figures, the physical parameters of these two ob- 
jects span a relatively broad range of values. In age and 
the two objects are quite distinct (Sll 5236 yielding older 



ages and a higher masses), reflected by the two separate 
peaks in the pdfs. For Sll 5236 the pdf is skewed toward 
low current SFR values, as the model favours short star- 
formation timescales and "old" ages. For these cases, how- 
ever, only weak Lya emission is expected, in contrast with 
the observations. Assuming a long star-formation timescale 
this shortcoming can be avoided; this situation is e.g. illus- 
trated by models with constant SFR (red dashed curves), in 
which case the SFR and extinction are found to be somewhat 
larger, and the stellar mass somewhat lower for Sll 5236. 
The fact that short star-formation timescales are favoured 
(for both objects, but more so for Sll 5236) also explains 
the relatively small overlap between the pdfs computed for 
variable r models (including the case of constant SFR) and 
the constant SFR case. 

For SGP 8884, the relatively high SFR and low mass 
implies a large specific SFR (sSFR> 100 Gyr~^) if taken at 
face value. In any case, relatively large uncertainties on the 
physical parameters remain, in particular given the small 
number of observational constraints. 

We now consider for comparison three other z — 5.7 
LAE for which more photometric constraints are available, 
and discuss the physical properties of these five LAEs to- 
gether. 



4.2 Comparison with SED modelling of other 
z = 5.7 LAE 



The three LAEs analysed earlier by Lai et al. (20071 are 



probably still the best comparison sample for our two 2 
bright LAEs. These objects were also selected by their 
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narrow-band excess, then spectroscopically confirmed, and 
finally selected for being detected in the IRAC 3.6 and 
4.5 /urn images. We have taken the published photometry 
from Lai et al. ( 2007 1 and we have reanalysed these objects 
using the same code and assumptions described abovsrl As 
an example, we show the best-fit SED of LAE #08 in Fig. 
[s] Clearly the effect of nebular emission (both lines and con- 
tinua) on the best-fit SED are readily visible. For all three 
objects, fits with nebular emission yield better fits (lower 
X^) than with standard templates. 



As already illustrated in Lai et al. (20071 the SEDs of 



their other objects are very similar, both in relative and ab- 
solute fluxes. Therefore it is not surprising to find a simpler 
probability distribution function for the physical parame- 
ters of Lai's objects (see blue dotted lines in Fig. |9| than 
the more bimodal pdf for SGP 8884 and Sll 5236. We find 
that the properties of SGP 8884 resemble more the typical 
values of the LAEs from Lai's sample, e.g. the SED, age, 
mass, SFR, and also Eb-y ■ 

Clearly do we find evidence for non-zero dust attenu- 
ation in all 5 high-z LAEs, as already found earlier by|Lai 



et al. (20071 for their objects and e.g. by Schaerer & Pello 
(20051 for other 2 > 6 LAEs. For example, for Lai's sample 
we find Eb-v ~ 0.15±0.05 and for our two objects -Bb-v ~ 
0.00-0.11 within 68% confidence. 

For the average stellar ages we do find a wide range 
of acceptable values (see Fig. [9|. Overall, relatively young 
ages (< 50-100 Myr) are however favoured, with a median 
slightly below 10 Myr for Lai's sample. Given these rela- 
tively young ages, the results from our SED fits are quite 
insensitive to the exact assumed SF history. 

The derived stellar masses also span a relatively large 
range, typically between few times 10* to few times 10^" 
Mq, the largest masses being found for Sll 5236. For Lai's 
sample we find a median stellar mass of ~ 6 x 10* M©, 
with 68% confidence limits b etween 2.2 and 6. 2 xlO* Mq. 
Compared to the analysis of |Lai et ah] ( |2007| ) our masses 
are significantly lower, by a factor 10 approximately. The 
difference is due to several effects in the SED fits: primarily 
the non-inclusion of nebular emission (which dominates the 
flux at 3.6 /im and explains a factor ~ 3-4 difference), and 
differences in the absolute scaling of the best-fit SEDs re- 
sulting from their peculiar treatment of the non-detections 
in the 5.8 and 8.0 /xm bands (~ factor 2.5). Our masses and 
other fit parameters should be more reliable as those from 
|Lai et alT| ( |2007[ ) since our models include the effects of nebu- 
lar emission, which must by selection be particularly strong 
for LAEs. 

The derived pdf for the SFR (Fig. [9) bottom right) 
sho ws aga in a rather large spread. For the objects from |Lai"et 



al. 



2007) we find SFR= UOtHl ^0 y^"^ (68% confidence 
interval), for our two objects a more complex distribution. 
The "standard" Lya (UV) calibrations!^ yield SFR ^ 23-30 
(20) Mq yr"^ for our two objects, without extinction correc- 
tions. Except for the case of variable star-formation history 
for Sll 5236 (cf. above) the (instantaneous) SFR derived 



^ We have noted that for these objects the pdf of the physical 
parameters does not change significantly if we allow for variable 
star-form ation histories or a ssu me constant SFR. 
^ See e.g. iKennicuttl (|1998| andlAjiki et al.|(|2003|l for Lyo. 



from the SED fits are higher due to non-zero extinction and 
due to ages younger than the timescale (~ 100 Myr) assumed 
for the UV calibration. 

Finally it is interesting to note the large speciflc star- 
formation rates (sSFR > 10-100 Gyr~^) found in general for 
the objects analysed here, except for the fits using variable 
star-formation history models for Sll 5236. Below we discuss 
the sSFR values in comparison with other data. 



4.3 Comparisons with other data from the 
literature 

The values of the extinction and stellar mass derived for 
the bulk of the LAEs (except maybe Sll 5236) agree well 
with those found from a larger sample of 2 « 6-8 LBGs by 



Schaerer fc de Barros|pOTo| ). The specific SFR (SFR/M*) of 
the LAEs appears, however, somewhat higher than average. 



Ono et al. (20101 have constructed the stacked SED of 
a large sample of LAEs at z — 5.7 and 6.5, which are unde- 
tected at 3.6 /im by IRAC. They have determined the aver- 
age physical properties from SED fits using models very sim- 
ilar to ours, including also nebular emission. For their sam- 
ple, showing a median Lya luminosity of (3.9 ± 0.27) x lO*'^ 
erg s~^, they find a low stellar mass (M*~ lO'^'^ M©), and 
no extinction. These differences with our estimates are most 
naturally explained if on average more massive galaxies suf- 
fer from more extinction, as well known from lower redshift 
(e.g., IBuat et al.|[2005l [2 008 Burgarcl la et al.|p007l [Daddil 
et al.|)2007| |Reddy et aLn2006, fiOOS) and also s uggested to 
hold at the highest redshifts from observations ( Schaerer & 
de Barros||2010[ ) and simulations ( |Dayal, Ferrara, fc Saro" 
20101. If Lya does not suffer much more strongly from ex- 



tinction, it may then also be expected, that the brightest 
Lya emitters at high-z are also among the most massive 
galaxies, at least on average. This is e.g. supported by the 



recent semi-analytical models of Garel et al. (20111, which 
use the results from Lya radiation transfer calculations to 
predict the Lya escape fraction depending on each galaxy's 
properties. 

We also find that our median -Bb-v value (_Eb-v ~ 



0.15 ± 0.05) determined for the LAE sample of Lai et al 
(2007) is in very good agreement with the average continuum 
attenuation fc = 0.23 predicted for LAEs at z = 5.7 by 



the models of Dayal, Hirashita, & Ferrara (20101. It is also 



basically identical to that adopted in the cosmological model 
of |Nagamine et aL] ( |2010] ). However, the above -Bb-v value 
is larger than a "cosmological average" of the continuum 



attenuation at 2 = 5.7 estimated by Hayes et al. (20111. 
Again, this is probably not surprising since our extinction 
value has been derived from relatively bright Lya emitters, 
which are also detected in the IRAC bands. How our median 
values thus relate to other averages is not trivial. 

In any case, our analysis of 5 of the brightest LAEs at 
2 ~ 5.7 confirms earlier flndings that some of these very 
redshift galaxies should contain some dust, and hence also 
be detectable in the infrared/sub-mm domain in the very 



near future ( Schaerer & Pell6|2005 


Boone et al.|2007 


Lai et 


al.|2007||Finkelstein et al.|2009||Dayal, Hirashita, & 1 


^errara 


2010 Dayal, Ferrara, & Saro|2010l. 



Finally, we note that the specific star-formation rate 
sSFR > 10-100 Gyr"^ obtained here for the LAEs is signif- 
icantly larger than the typical value of ~ 2-3 Gyr~^ found 
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for LBGs by other authors at z ~ 2 and at redshifts up 



SGP 8864 



to .z > 6 (see e.g. Daddi et al. 2007 Stark et al. 2009 



Gonzalez et al. 2010 McLure et al. 


20111. However, some 


LBGs at 2 ~ 2-3 (Erb et al.|2006b|c 


Yoshikawa et al.|2010 



Dessauges-Zavadsky et al.||2011 i, some LAEs at z ~ 3 — 4 



(Ono et al. 2010L and some LBGs at z > 4 (Schaerer & 



de Barros 2010 de Barros et al. 20111 also show compa- 



rably high sSFR (and a large scatter) Compared to recent 
semianalytic and hydrodynamic galaxy formation models a 
low sSFR ~ 2-3 Gyr~^ constant with redshift above z > 2 



seems problematic (cf. 


Drory & Alvarez|2008, 


Nagamine et 


al.||2010| IBouche et al.||2010| |Khochfar et al. 


|2011| IWein- 



appear larger than predicted by these models, although they 
generally predict an increase of sSFR with decreasing galaxy 
mass, which might qualitatively explain why the LAEs anal- 
yses here have high sSFR. More in-depth SED studies of well 
constrained LAEs and LBGs at different redshift and inde- 
pendent measurements of SFR and mass will be needed to 
clarify these issues further. 




elength [^tm] 



4.4 Dust extinction 

The modelling of the SED and the modelling of the Lya line 
provide independent estimates of the amount of absorption 
from dust. While the allowed values from the modelling of 
the Lya line are relatively unconstrained, the estimates of 
the amount of absorption derived in this way are generally 
different to the estimates derived from the modelling of the 
SED. In particular, the estimate of the amount of dust ab- 
sorption in SGP 8884 is a factor of 2 to 3 higher than that 
obtained from the SED fit. The difference has been noted 
in other LAEs ( [Verhamme et al.|2008"l [Vanzella et al.|2010| ) 
and may be related to the difference that has been seen lo- 
cal star- forming galaxies ( Calzetti|2001 1 , since one measure 
is based on the fit to Lya, which is entirely nebular emis- 
sion, and the other is based on the fit to the SED, which is 
mixture of stellar light, nebular emission and nebular con- 
tinuum. Alternatively, this apparent difference may also be 
related to the simplified model assumptions in the radiation 
transfer code (e.g. simplified geometry) or in the SED fits. 

The Lya fits to both SGP 8884 and Sll 5236 also result 
in unusually low gas-to-extinction ratios, with values around 
Nu/E{B - V) ~ 3 - 6 X 10^^ cm~^ mag"^ for the models 
that exclude the effects of the IGM. Compared to the mean 
ratio in the Galaxy, Nh/E{B-V ) = 5.8 x 10^^ cm"^ mag"\ 
( Bohlin, Savage, fc Drake] 19781 the gas-to-extinction ratios 
in SGP 8884 and Sll 5236 are about 100 times lower. The 
values are also lower than those estimated for other LAEs 
( Verhamme et aL]|2008 Vanzella et al.pOlO l. However, the 
range allowed by the fits is also large. For example, the 68.3% 
confidence interval on the extinction allows values as high 
Nii/E{B-V) = 4x 10^° cm-2 mag"^ for SGP 8884, which is 
still lower than that measured in the Galaxy, but similar to 
the values estimated for other LAEs. As already mentioned, 
it remains open whether these results can be taken at face 
value or whether they are due to simplifying assumptions in 
our model of modelling of the Lya line. 



Figure 6. Observed (blue symbols) and best-fit SED of SGP 8884 
using templates including nebular emission (black lines and red 
symbols showing the synthetic fluxes) and standard spectral tem- 
plates (magenta dashed line) for solar metallicity. For the other 
metallicities, the fits are very similar. The upper limits (1 a val- 
ues) are also included in the SED fits. 



5 SUMMARY 

We obtained deep optical and near-IR spectra and Spitzer 
IR imaging of two of the most luminous LAEs known at 
z = 5.7 and we use these data to constrain the properties of 
these galaxies. 

The continuum red-ward of Lya is clearly detected in 
both objects, allowing us to derive reasonably precise es- 
timates of the rest-frame equivalent widths. For both ob- 
jects we find equivalent widths that are around 160 A with 
uncertainties of 10 and 20% for SGP 8884 and Sll 5236, 
respectively. 

We used an analytic model and a 3-D Monte Carlo ra- 
diation transfer code with two treatments for how the IGM 
affects the line to model the line profile. All models, accu- 
rately model the blue edge of the line and the continuum 
out to velocities of > 1000 kms~^ blue-ward of Lya. At a 
resolution of 5000 and with the present S/N, we cannot dis- 
tinguish the three models. Consequently, from the blue side 
of the Lya line we do not find observational proof that Lya 
of our two LAE is affected by the IGM at these redshifts, al- 
though the continuum is found to be depressed over a larger 
wavelength range blue- ward of Lya. 

However, neither the analytic model nor the model that 
excludes the IGM from affecting the line profile adequately 
model the red wing of the line in either LAE. Only the ra- 
diation transfer model that includes a simple prescription 
for how the IGM affects the line adequately models the red 
wing in both LAEs. While this may be used as evidence to 
support the idea that Lya is being affected by the IGM at 
these redshifts, it is not yet conclusive, as the distribution 
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Figure 9. Probability distribution functions (pdf) for the derived physical parameters of the two LAEs SGP 8884 (black) and Sll 5236 
(red lines) and the 3 2 = 5.7 LAEs from [Lai et al.| | |2007| l (blue dotted lines). For our objects two pdfs are shown, assuming variable 
star-formation histories (solid lines) or constant SFR (dashed); the former underpredicts the observed Lya equivalent width of Sll 5236. 
For Lai's objects the pdf considers variable star-formation histories. The area under the pdfs are normalised to a common value for the 
Lai et al. sample and for our objects taken together. Note that for Sll 5236 the pdf also extends to lower SFR than shown on the figure. 
Top left: Attenuation Eb— v derived for the |Calzetti| j2001[ | law. Top right: Stellar age. Middle left: stellar mass. Middle right: 
instantaneous star-formation rate. Bottom: specific star-formation rate. 
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Figure 8. Observed (blue symbols) and best-fit SED of object 
#08 from |LaietaI] ( |2007| l using templates including nebular emis- 
sion (black lines and red symbols showing the synthetic fluxes) 
and standard spectral templates (magenta dashed line). The up- 
per limits (1 cr values) are also included in the SED fits. 



Figure 7. Observed and best-fit SED of Sll 5236 using the same 
symbols as in Fig. [6] A constant SFR has been imposed for the 
fit with nebular lines to assure a sufficiently strong Lya emission 
(cf. text). 



of dust and neutral hydrogen in the ISM in these LAEs is 
likely to be more complex than assumed in the models used 
here. 

The radiation transfer models predict that the intrinsic 
equivalent width is about double the observed one, about 
300 A which is at the upper end of the range allowed for 
a young, moderately metal-poor star-forming galaxy. This 
is independent of how we have treated the IGM. However, 
uncertainties are currently large and values as high as 700 A, 
which in the realm expected for a Population III burst, are 
allowed. 

Both LAEs were observed and detected with IRAC on 
the Spitzer Space Telescope. We combined the IRAC pho- 
tometry with a measurement of the continuum at 1300 A 
and ground based data to constrain the SED of these two 
objects and compared them to 3 other LAEs at 2 = 5.7 from 
the literature ( |Lai et al.||2007[ ). The SEDs of these 5 LAEs 
have been analysed with state-of-the-art SED models includ- 



ing the effects of nebular emission (cf. [Schaerer &: de Barro8| 



2009 2010 1. In terms of age and mass, our two LAEs are 



quite distinct from each-other. Sll 5236 appears older and 
more massive than SGP 8884. Four of the five LAEs anal- 
ysed here show masses of the order of ~ 5 x 10* M0, lower 
than previous estimates of Lai et al. We find evidence for 
the presence of some dust in all objects, and indications for 
fairly high specific star-formation rates (> 10-100 Gyr~^). 



Nebular line emission, principally from [OIII] AA 
4959,5007 doublet and H/3, makes a significant contribution 
to the flux in th e IRAC 3.6 ^m fi lter. For SGP 8884 and the 
three objects in |Lai et al.| ( [20071 ), dominates. Dusty, and 
rich in emission lines, these objects will be prime targets 
for the next generation of extremely large telescopes, JWST 
and ALMA. 
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